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The neutral B meson lifetime is measured with the data collected by the BABAR detector at the 
PEP-II storage ring during the years 1999 and 2000, with a total integrated luminosity of 20.7 fb~ . 
The decays B° -> D*~tt+ and B° -» D*~p + are selected with a partial-reconstruction technique, 
yielding samples of 6970 ± 240 and 5520 ± 250 signal events, respectively. With these events, the 
B° lifetime is measured to be 1.533 ± 0.034 (stat.) ± 0.038 (syst.) ps. This measurement serves as a 
test and validation of procedures required to measure the CP violation parameter sin(2/3 + 7) with 
partial reconstruction of these modes. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The neutral B meson decay modes jlj B° — > D*~h + , where h + is a light hadron (n + , p + , a^), have been pro- 
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posed for use in theoretically clean measurements of 
sin(2/3 + 7) Q , where (2/3 + 7) is a combination of angles 
of the Cabibbo-Kobayashi-Maskawa fl unitarity trian- 
gle. Since the time-dependent CP asymmetries in these 
modes are expected to be of order 2%, large data samples 
and multiple decay channels are required for a statisti- 
cally significant measurement. The technique of partial 
reconstruction of D*~ mesons, in which only the soft pion 
ir s from the decay D*~ — > D°ir~ is reconstructed, has al- 
ready been used to select large samples of B meson can- 
didates 0. This technique is applied here to the decays 
B° —> D*~ir + and B° — > D*~ p + in order to measure the 
B° lifetime. This analysis constitutes a first step toward 
measuring sin(2/3 + 7), validating the procedures devel- 
oped for candidate reconstruction, background character- 
ization, vertex reconstruction, and fitting of decay time 
distributions. These procedures address the main com- 
plications introduced by partial reconstruction, namely 
the large background and the tracks originating from the 
unreconstructed D°, which may affect the vertex recon- 
struction. 

The analyses applied to the B° — » D*~tt + and 
B° — > D*~ p + modes are similar. Detailed differences 
between them are the result of optimization in the pres- 
ence of the different background characteristics in the two 
modes. Additional details regarding the analysis proce- 
dures can be found in Refs. || and ||. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy stor- 
age ring during the years 1999 and 2000. The data con- 
sist of 22.7 million BB pairs, corresponding to an in- 
tegrated luminosity of 20.7 fb _1 recorded at the T(AS) 
resonance. In addition, 2.6 fb" 1 of "off- resonance" data 
were collected about 40 MeV below the resonance. Sam- 
ples of simulated BB and continuum e + e _ — > qq events, 
where q stands for a u, d, s, or c quark, were gener- 
ated using a GEANT3-based detector simulation Q and 
processed through the same reconstruction and analy- 
sis chain as the data. The equivalent luminosity of the 
simulated events is approximately one third the data lu- 
minosity. We also used signal Monte Carlo samples with 
an equivalent luminosity several times larger than that 
of the data. 

The BABAR detector, described in detail elsewhere [||, 
consists of five subdetectors. Charged particle trajecto- 
ries are measured by a combination of a five-layer sil- 
icon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH) in a 1.5 T solenoidal magnetic field. Tracks 
with low transverse momentum are reconstructed by the 
SVT alone, thus extending the charged particle detec- 
tion down to transverse momenta of ~ 50MeV/c. Pho- 
tons and electrons are detected in a CsI(Tl) electromag- 
netic calorimeter (EMC), with photon energy resolution 
a E /E = 0.023( J B/GeV)- 1 / 4 © 0.019. A ring-imaging 
Cherenkov detector (DIRC) is used for charged parti- 
cle identification. The instrumented flux return (IFR) is 



equipped with resistive plate chambers to identify muons. 

In the partial reconstruction of a B° — > D*~h + can- 
didate, only the hadron h and the ?t s tracks are recon- 
structed. The angle between the momenta of the B and 
the h in the center-of-mass (CM) frame is then computed: 

M 2 D , _ - M% o - Ml + E CM E h 
cos B Bh = —2 2 — * 1 

2PB\Ph\ 

where M x is the mass of particle x, Eh and ph are the 
measured CM energy and momentum of the hadron h, 
Ecm is the total CM energy of the beams, and pb = 

\J ^cm ~ M b° • ^ u masses refer to the nominal val- 
ues B, except in the case h = p, where the measured 
7r + 7T invariant mass m(7r + 7r°) is used. Events are re- 
quired to be in the physical region | cos 9 Bh\ < 1- Given 
cos 8 Bh and the measured four- momentum of h, the B 
four-momentum can be calculated up to an unknown 
azimuthal angle </> around ph- For every value of 0, 
the expected D° four-momentum Vd{4>) is determined 
from four-momentum conservation, and the (^-dependent 
"missing mass" is calculated, m(<j>) = ^/\Vd{4>)\ 2 ■ We 
define the missing mass m m j SS = = [m max + m m in] , where 
m max and m m i n are the maximum and minimum values 
of m(</>). In signal events, this variable peaks at the nom- 
inal D° mass M D o, with a spread of about 3 MeV/c 2 for 
B° -> D*-n+ (3.5 MeV/c 2 for B° -> D*~p+) filj, while 
the distribution of background events is broader. The 
missing mass is the main variable used to distinguish sig- 
nal from background. 

We define the D* helicity angle 9d* to be the angle 
between the directions of the D° and the B° in the D* 
rest frame. This variable is used in the event selection 
described below. In the B° — > D*~n + analysis, 6jj, 
is computed assuming that the B momentum lies in the 
plane defined by the h and ir s momenta in the CM frame. 
This assumption also yields the D° direction. In the 
B° — » D*~ p + analysis, the value of cos#d» is computed 
by applying the constraint m m i ss = M D o giving two pos- 
sible solutions for the D° direction M. In B° — > D*~ p + , 
the p helicity angle 9 P is defined as the angle between the 
directions of the ir° (from the decay of the p) and the 
CM system in the p rest frame. 

We select events in which the ratio of the 2nd to the 
0th Fox- Wolfram moment [pj, computed using charged 
particles, is smaller than 0.35. The candidate B° daugh- 
ter tracks are required to originate within 1 cm (1.5 cm) 
of the interaction point in the x-y plane (the plane per- 
pendicular to the beams), and within ±4 cm (±10 cm) 
of the interaction point along the direction of the beams. 
Tracks are rejected if they are highly likely to be a kaon 
or a lepton on the basis of their ionization, Cherenkov 
angle, energy deposited in the EMC, and pattern of hits 
in the IFR. 

B° — > D*~n + candidates are rejected if another track 
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is found within 0.4 rad of the momentum of the hard pion 
7T/j [Q in the CM frame. This requirement helps to reject 
continuum events, where tracks tend to be clustered in 
jets. A Fisher discriminant |l3| is computed from 
15 event shape variables. Among these variables is the 
scalar sum of the CM momenta of all tracks and neutral 
candidates in nine 20° single-sided cones around the Tth 
direction. We require |cos6*d*| to be larger than 0.4. A 
cut on Fn is used to reduce the continuum background. 

In the reconstruction of B° — ► D*~ p + candidates, the 
charged p candidates are identified by their decay to 
a hard charged pion ir^ and an . To suppress fake 
7T° candidates, the 7r° momentum in the CM frame is 
required to be greater than 400MeV/c. The invari- 
ant mass of the 7r° — > 77 candidate must be within 
20 MeV/ c 2 of the nominal ir° mass |9| . The invariant mass 
?7i(7T + 7r ) of the p candidate must be between 0.45 and 
1.10 GeV/c 2 . To suppress combinatoric background, we 
require | cos6> p | > 0.3 and | cos 9d* \ > 0.3, and also reject 
events that satisfy both cos6> p > 0.3 and cos6*£>* < —0.3. 
A Fisher discriminant F p is computed using the scalar 
sum of the CM momenta of all tracks and neutrals in nine 
10° double-sided cones around the p direction. In about 
10% of the events, more than one partially reconstructed 
candidate per event satisfies all the requirements. In 
such events only the candidate with the smallest value 
of |m m i SS — M D o \ in the event is used. 

The decay position z roc of the partially reconstructed 
B candidate along the beam direction is determined by 
constraining the tt^ and the tt s tracks (only the Wh track 
for B° — ► D*~ p + ) to originate from the beam-spot in 
the x-y plane. The beam spot is determined on a run- 
by-run basis using two-prong events Its size in the 
horizontal direction is 120 /im. Although the beam spot 
size in the vertical direction is only a few microns, a beam 
spot constraint of 30 /jm is applied, so as to account for 
the flight of the B° in the vertical direction. 

The decay position z ot her of the other B meson along 
the beam direction is measured with all tracks, exclud- 
ing Wh, 7r S j and any track whose CM angle with respect 
to the D° direction (either of the two calculated direc- 
tions in the B° — ► D*~ p + case) is smaller than 1 radian. 
This "cone cut" reduces significantly the number of D° 
daughter tracks used in the other B vertex. The tracks 
satisfying this requirement are fit with a constraint to the 
beam-spot in the x — y plane. The track with the largest 
contribution to the % 2 of the vertex, if greater than 6, is 
removed from the vertex, and the fit is carried out again, 
until no track fails this requirement. B° — » D*~ir + can- 
didates are required to have at least two tracks remaining 
in the other B vertex. 

The z distance between the two B decay vertices, 
Az = Ztec — Mother, is computed. Fitting the resid- 
ual Az — Aztruc in simulated events, where Aztrue is the 
true Az, with the sum of two Gaussians, we find that 
67% (57%) of the B° -> D*~7r+ (B° -> D*~p+) events 



lie in the core Gaussian of width 116 /im (178 pm). The 
Az; resolution is dominated by the measurement of z therj 
and by the z rcc measurement when the ith transverse mo- 
mentum is below about 400 MeV/c. 

The decay time difference At is then calculated using 
the approximation At « Az/ (7/fc), where the CM frame 
boost 7/3 is determined from the beam energies, and has 
an average value of 0.55. This approximation results in 
a 0.2 ps r.m.s. spread in the calculation of At. 

For B° — » D*~ir + candidates, At is computed apply- 
ing an event-by-event correction to the measured value 
of Az. This correction, determined from the simulated 
signal sample as a function of Az, removes the bias in 
Mother due to the tracks coming from the D° decay. With- 
out correction, the effect of this bias would be to re- 
duce the measured lifetime by approximately 4%. In the 
B° — > D*~ p + analysis a different correction is applied to 
the measured lifetime value, as explained later. 

The estimated error <7 At in the measurement of At is 
calculated from the uncertainties in the parameters of 
the tracks used in the two vertex fits. A requirement on 
the vertex fit probabilities removes badly reconstructed 
vertices. For both modes we also require Ai| < 15 ps 
and er At < 2.4 ps (er At < 4 ps for B° — > D*~p + ). 

After applying all the above requirements, we find 
four broadly-defined types of events that contribute to 
the background: (1) continuum events; (2) combina- 
toric BB background due to random h and 7r s combi- 
nations; (3) B° -> D*~p+ (B° -»• D*~a+) decays in the 
B° -> D*~tt + (B° -> D*~p+) sample; (4) peaking BB 
events, which are distributed as a broad peak in the m m j SS 
spectrum. The peaking background is mostly due to 
B — ► D**n decays in the B° — > Z)*~7r + sample. In the 
B° — > D*~ p + sample, it is due to signal events in which 
the 7T/J candidate originates from the other B. 

The lifetime t b o is obtained from an unbinned maxi- 
mum likelihood fit, as described below, with a probability 
density function (PDF) J- (At, <r At , £). Here £ refers to 
the kinematic variables used to distinguish signal from 
background. For B° — > D*~ir + we set £ = m m ; ss ; for 
B a ->D*-p+ we set £ = (m miss , m(7r + 7r°), F p ). The 
PDF has the form 

?((,, At, (T Af ) = /cont^Ccont (O-^cont (At, (T Af ) 

~t~ /comb^comb (O^comb (At, a At ) 
+ fD*xKD-x(QF D >x(At,<j^) 

+ /poak^Cpcak(£)^ r poak(At, (T At ) 

+ / sig /C sig (£)^ sig (At,a At ), (2) 

where subscripts refer to the four types of backgrounds 
enumerated above and to signal events. For each event 
type i, ft is the relative population of these events in the 
data sample, JCi(£) is their kinematic- variables PDF, and 
Ti(At, <7 At ) is their time-dependent PDF. The constraint 
ft = 1 is enforced. 
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For B° — > D*~ir + , /Q(m m i ss ) consists of binned histo- 
grams obtained from the Monte Carlo simulation. For 
B° — > D*~p + candidates, we use the product = 
M l (m m{ss )1l i {m('K + TT Q ))'D i (F p ) 1 where M l {m miss ) is the 
sum of a bifurcated Gaussian and an ARGUS func- 
tion Q, 7^i(m(7r + 7r )) is the sum of a parabolic back- 
ground and a relativistic P-wave Breit-Wigner, and 
T>i(F p ) is a bifurcated Gaussian function. 

For each event type i, !Fi(At, a At ) is the convolution 
TV / P(At tIUC )R({At - At u . uc )/a At )dAt tIUC of the "true" 
distribution P(Ai tr uc) and the detector resolution func- 
tion R((At — Aitrue)/c At ), which is parameterized as 
the sum of three Gaussian distributions. N is a nor- 
malization constant. The parameters of P(Ai t ruc) and 
R((At — AttrueV^A*) are obtained separately for each 
event type. For signal events of both modes we take 
P(A^truo) = e l Attiuc l/ 7 ^o This functional form is also 
used for the combinatoric and peaking backgrounds in 
B° — » £>*~7r + , but with independent parameters. In 
B° — > D*~ p + , the source of the peaking background mo- 
tivates its distribution to be P(At tr uc) = ^(Attruo), and 
the distribution used for the combinatoric background is 
P(Ai truo ) = ae -l A **.-l/^' + (i _ a )S(At tIue ), with an ef- 
fective lifetime parameter r'. J-D*x(At, <x At ) is assumed 
to be identical to J- s i g (At, <J At ). The continuum back- 
ground is modelled as P(At tr uc) = &e~l Attel »l/' ro ° nt + (1 - 
b)d{At tmc ). 

Several subsamples are defined and used in the life- 
time fit. Events with a candidate in which the h and tt s 
have opposite charges and with m miss > 1.860 GeV/c 2 
(rn miss > 1.845 GeV/c 2 in B° D*~ p+) constitute 
the "signal region" sample. Those satisfying 1.820 < 
m miss < 1.850 GeV/c 2 (1.810 < m miss < 1.840 GeV/c 2 ) 
constitute the "sideband" . Events in which h and tt s 
have the same charge are labeled as "same-charge". In 
the B° — » D*~ix + analysis, we apply a requirement on 
the Fisher discriminant that suppresses BB events, to 
select a "PP-depleted" sample that is enriched in con- 
tinuum events. The sideband, same-charge, and BB- 
depleted samples serve as control samples for studying 
the At distributions of the backgrounds. 

In the B° — > D*~ p + analysis, about 11.5% of the par- 
tially reconstructed signal events are also fully recon- 
structed in the D° decay modes D° — > K + tt~ or 
K + ir~ir°. This yields a sample that, while relatively 
small, has a low background contamination of about 5%. 
This clean signal sample is used in the fits described be- 
low, improving the determination of the signal PDF pa- 
rameters. 

The B° lifetime t a is obtained in a three-step proce- 
dure using signal region and control sample events. 

In the first step, the fractions fi in the signal re- 
gion and in the different control samples are obtained 
from kinematic- variable fits conducted simultaneously on 
the on- and off-resonance samples (and the fully recon- 
structed sample for the B° — > D*~ p + signal region). The 



fit PDF is that of Eq. (0), but with all ^(Af,cr At ) re- 
placed by unity. In the B — > D*~n + analysis this fit de- 
termines /peak and /cont- The fraction of B° — » D*~ p + 
events fwx in the B° — ► D*~tt + sample is assumed to 
be 16.8%, as predicted by the Monte Carlo simulation 
and the relative branching ratio || . This fit (Fig. [l](a)) 
yields 6970 ± 240 signal B° -> D*~7r+ events. In the 
B° — > D*~ p + analysis the kinematic-variable fit deter- 
mines /cont, as well as all the parameters of /C C ont(£), 
A / isig(w m i ss ), and 7^. s i g (r7i(7r + 7r )). The parameters of 
£>sig(P P ), £ comb (£), and /C pea k(£), as wel1 as / P oak// s i g 
(9.7%) and fr>* ai /fsig (11.6%), are obtained from the 
Monte Carlo simulation. The kinematic-variable fit to 
the B° -> D*~p+ sample (Figs. |(a), (b) and (c)) yields 
5520 ± 250 B° -> D*~p + events, including 691 ± 36 fully 
reconstructed events. 

In the second step, all the parameters determined in 
the first step are fixed, and the parameters of J-i(At, a At ) 
of the backgrounds are determined entirely from the con- 
trol data samples. In the B° — > D*~ir + case, the pa- 
rameters of ^cont(At, u At ) are obtained from a fit to the 
PP-dcplcted sample, and those of the J 7 C omb(At, a At ) 
are obtained from the same-charge sample. The pa- 
rameters of ^-p 0a k(At, <7 A1 ) are assumed to be identical 
to JT comb (At, cr At ). In P° — > D*~p + , the parameters of 
•? r comb(Ai, (T At ) are determined from the sideband sam- 
ple, and those of J-" pea k(Af, <r At ) are obtained from the 
same-charge sample. Each of the B° — > D*~ p + control 
sample fits is conducted simultaneously on the on- and 
off- resonance data, and the parameters of T con t{At, cr At ) 
are determined for each control sample simultaneously 
with the BB PDF parameters. 

In the final step, using the background J-^At, cr At ) pa- 
rameters obtained in the previous step, the signal region 
sample is fit to extract the signal J- s i g (At, a At ) param- 
eters. In B° — ► D*~n + this fit has six free parameters 
describing J r s i g (At, <r At ). In P° — > D*~ p + , the fit is done 
simultaneously to on- and off-resonance events, as well 
as fully reconstructed events, and has 15 free parameters 
describing _F sig (At, cr At ) and J" cont (At, a At ). 

The results of the last fit step, shown in Figs. |l|(b) 
and|(d), are t = 1.510 ± 0.040 ps for B° -> D*~it+ 
and t b0 = 1.616 ± 0.064 ps for B° -> D*~p + , where the 
errors are statistical only. These results are obtained af- 
ter a correction of -0.014±0.020 ps (+0.071 ±0.028 ps for 
B° —> D*~ p + ), determined from the Monte Carlo simu- 
lation. The correction accounts for biases due to the fit 
procedure, the event selection and, in the B° — > D*~ p + 
case, the effect of D° daughter tracks passing the cone 
cut and being used for the determination of the other B 
vertex. The errors in the corrections are propagated to 
the final result as systematic errors. 

The systematic uncertainties are listed in Table |, and 
described here. (1) The fractions and the PDF pa- 
rameters of the background components were varied by 
their statistical errors, taking into account mutual cor- 
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FIG. 1: Distributions of (a) missing mass and (b) At for can- 
didate B° — > D*~tt + events. The result of the fit (solid line) 
is superimposed on data (data points). The hatched, cross- 
hatched and shaded areas are the peaking BB, combinatoric 
BB, and continuum contributions, respectively. The At plot 
is obtained with the requirement m m i ss > 1.860 GeV/c 2 . 



relations, obtained from the fits of the first two analy- 
sis steps. (2) The PDF parameters and lifetime correc- 
tions that were obtained from the Monte Carlo simu- 
lation were varied by the statistical error in the Monte 
Carlo fits. The full analysis chain, including event re- 
construction and selection, was tested with the Monte 
Carlo simulation, and the statistical precision of the con- 
sistency between the generated and fit lifetimes was as- 
signed as a systematic error. The Monte Carlo statis- 
tical errors in the evaluation of the various corrections 
described above were propagated to the final result. (3) 
The level of B° -> D*~ p + (B° -> D*~a% ) background in 
the B° — > D*~tt + (B° — > D*~ p + ) sample was varied by 
the relevant branching fraction errors Q , and the fraction 
of B — > D**p + background events in the B° — > D*~ p + 
sample was varied between and 40% of the signal yield. 
(4) The number of D° tracks satisfying the cone cut in 
the simulated sample was varied by ±5% and the asso- 
ciated bias was reevaluated. (5) The parameters of J- S i g 
that were fixed in the fits were varied within conservative 
ranges. (6) Extensive parameterized Monte Carlo simula- 
tion studies were conducted to evaluate statistical biases 
in the fits due to limited data sample size or as the result 
of changes in the functional form of R((At— Ait me )/<7 At ). 
(7) The At fit range was varied between \At\ < 10 ps 
and | At | < 20 ps. (8) The z length scale of the de- 
tector has been determined with an uncertainty of 0.4% 
from the reconstruction of secondary interactions with a 
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FIG. 2: Distributions of (a) missing mass, (b) p candi- 
date invariant mass, (c) Fisher discriminant F p and (d) At 
of B° — » D*~ p + candidate events. The result of the fit (solid 
line) is superimposed on data (data points). The hatched, 
cross-hatched and shaded areas are the peaking BB, combina- 
toric BB, and continuum contributions, respectively. The At 
plot is obtained with the requirement m m i ss > 1.854 GeV/c 2 , 
0.60 < m(7T + 7T°) < 0.93GeV/c 2 , and F p < -2.1. 



TABLE I: Summary of the systematic uncertainties on the 
measured B° lifetime. 



Source 


Errors (ps) 




B °^D*-7v+ B°^D*-p+ 



(1) Background parameters 


0.023 


0.044 


(2) Monte Carlo statistics 


0.021 


0.042 


(3) Fractional composition 


0.008 


0.024 


(4) D° tracks bias 


0.017 


0.026 


(5) At resolution model 


0.011 


0.015 


(6) Likelihood fit bias 


0.005 


0.016 


(7) At range 


0.009 


0.009 


(8) 2 scale 


0.006 


0.007 


(9) SVT misalignment 


0.008 


0.008 


(10) Beam energies 


0.002 


0.002 


Total 


0.041 


0.075 



beam pipe section of known length [|15[ . The systematic 
uncertainties related to the detector alignment (9) and 
beam energy uncertainty || (10) were also taken into ac- 
count. The total systematic error in the B° — > D*~ir + 
(B° -> D*~p+) analysis is 0.041 ps (0.075 ps). 

Several cross-checks were conducted to ensure the va- 
lidity of the result. The data were fit in bins of the lab 
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frame polar angle, azimuthal angle, and momentum of 
the TVh, and in subsamples corresponding to different SVT 
alignment calibrations. The fit was repeated with dif- 
ferent values of the cone cut ranging from 0.75 to 2.00 
radians (0.6 to 1.2 radians for B° — > D*~ p + ). Different 
functional forms of R((At ~ At truc ) / a &t ) were used in the 
fit . In all cases, no statistically significant variation of the 
result was observed, beyond those already accounted for 
in the systematic errors. 

In summary, in a sample of 22.7 million BB pairs, 
we identify 6970 ± 240 B° -> D*~n+ and 5520 ± 250 
B° — > D*~ p + partially reconstructed decays. These 
events are used to measure the B° lifetime, obtaining 



t 



1.510 ± 0.040 (stat.) ± 0.041 (syst.) ps in F° 



D*-n+ and t = 1.616 ±0.064 (stat.) ± 0.075 (syst.) ps 
in B° — > D*~ p + . The combined measurement, taking 
intoi account correlated errors, is 



r = 1.533 ±0.034 (stat.) ± 0.038(syst.) ps. 

This result is in good agreement with the world average 
B° lifetime r = 1.542 ± 0.016 ps § and with other 
recent BABAR measurements fl6f , confirming the validity 
of using partially reconstructed events in time dependent 
measurements. 
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